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Abstract 

We investigate cold dense matter in the context of Sakai and Sugimoto's holographic 
model of QCD in the Wigner-Seitz approximation. In bulk, baryons are treated as 
instantons on x in each Wigner-Seitz cell. In holographic QCD, Skyrmions 
are instanton holonomies along the conformal direction. The high density phase is 
identified with a crystal of holographic Skyrmions with restored chiral symmetry at 
about 4M^^/tt^. As the average density goes up, it approaches to uniform distribution 
while the chiral condensate approaches to p-wave over a cell. The chiral symmetry is 
effectively restored in long wavelength limit since the chiral order parameter is averaged 

5/3 

to be zero over a cell. The energy density in dense medium varies as , which is the 
expected power for non-relativistic fermion. This shows that the Pauli exclusion effect 
in boundary is encoded in the Coulomb repulsion in the bulk. 



1 Introduction 



The AdS/CFT approach [1] provides a powerful framework for discussing large Nc gauge 
theories at strong couphng A = g^Nc- The model suggested by Sakai and Sugimoto (SS) [2] 
offers a specific holographic realization that includes Nf flavors and is chiral. For Nf <C 
chiral QCD is obtained as a gravity dual to Nf D8-D8 branes embedded into a D4 background 
in 10 dimensions where supersymmetry is broken by the Kaluza-Klein (KK) mechanism. The 
SS model yields a holographic description of pions, vectors, axials and baryons that is in good 
agreement with experiment [2, 3, 4]. The SS model at finite temperature has been discussed 
in [5] and at finite baryon density in [6, 7, 8]. The finite baryon chemical potential problem 
in D3/D7 model was discussed in [9] and Isospin chemical potential and glueball decay also 
has been discussed in [10] . 

Cold and dense hadronic matter in QCD is difficult to track from first principles in current 
lattice simulations owing to the sign problem. In large Nc QCD baryons are solitons and a 
dense matter description using Skyrme's chiral model [11] was originally suggested by Skyrme 
and others [12]. At large Nc and high density matter consisting of solitons crystallizes, as 
the ratio of potential to kinetic energy T = V/K ^ N^ is much larger than 1. QCD matter 
at large Nc was recently revisited in [14]. 

The many-soliton problem can be simplified in the crystal limit by first considering all 
solitons to be the same and second by reducing the crystal to a single cell with boundary 
conditions much like the Wigner-Seitz approximation in the theory of solids. A natural 
way to describe the crystal topology is through with periodic boundary conditions. In 
so far, this problem can only be addressed numerically. A much simpler and analytically 
tractable approximation consists of treating each Wigner-Seitz cell as S''' with no boundary 
condition involved. The result is dense Skyrmion matter on [13]. Interestingly enough, 
the energetics of this phase is only few percent above the energetics of a more involved 
numerical analysis based on T^. Skyrmions on restore chiral symmetry on the average 
above a critical density. While Skyrmions on are unstable against T^, they still capture 
the essentials of dense matter and chiral restoration in an analytically tractable framework. 

Cold dense matter in holographic QCD is a crystal of instantons with F = VX/vp ^ 1 
where vp ^ 1/Ncis the Fermi velocity. (In contrast hot holographic QCD has F = y/X » 1). 
When the wigner-Seitz cell is approximated by S^, the pertinent instanton is defined on 

X R. In this paper, we investigate cold QCD matter using instantons on x R in bulk. 
As a result the initial D4 background is deformed to accomodate for the which is just the 
back reaction of the flavour crystal structure on the pure gauge theory. Holographic dense 
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matter can be organized in 1/A at large Nc- In our model the baryon density is uniform while 
the chiral condensate is p-wave over a cell. However, the chiral condensation is averaged to 
be zero over a cell so that the chiral symmetry is effectively restored in long wavelength limit. 
We will show that as the average density goes up, it approaches to uniform distribution while 
the chiral condensate approaches to p-wave over a cell. The energy density in dense medium 
varies as n^^^, which is the expected power for non-relativistic fermion. This shows that the 
Pauli exclusion effect in boundary is encoded in the Coulomb repulsion in the bulk. 

In section 2, we define this deformation and discuss the D8 brane embedding structure. 
The instantons on S*^ x i? in the flavour D8 brane is discussed in section 3. In section 4,5,6 
we derive the equation of state of cold holographic matter using the small size instanton 
expansion and in general. In section 7 we show how the holographic small instantons in bulk 
transmute large size Skyrmions on the boundary. The comparison to other models of nuclear 
matter is carried in section 8. Our conclusions are in section 9. 



2 D8 brane action 

We consider crystallized skyrmions at finite density in the Wigner Seitz approximation. 
Spatial is naturally converted to with periodic boundary conditions. As a result the 
D4 background geometry is deformed. The baryons are then instantons on x R. Most 
solutions are only known numerically on the lattice. A simpler and analytically tractable 
analysis that captures the essentials of dense matter is to substitute by in bulk with 
no boundary conditions altogether. As a result, the D4 background dual to the crystal is 
modified with the boundary special space as S^. Specifically, the 10 dimensional space is 
that of {R^ X S^) X R^ X S^. The ensuing metric on D4 is therefore 

ds' = (^) i-dt' + oi'dnl + mdr') + (-^ + u'dnl^ , (1) 

dQs = di/j^ + sin^ij d9^ + sm^jjj sin^O , f{U) = l- , (2) 

^'-dsi^^y, F4^dC3=^64, (3) 

While this compactified metric is not an exact solution to the general relativity (GR) equa- 
tions for small size S^, it can be regarded as an approximate solution for large size S^. 
Indeed, in this case, the GR equations are seen to be sourced by terms wich are down by the 
size of S^. Here, (3) can be regarded as an approximation to the stable metric with for 
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a dense matter analysis. Clearly, the former is unstable against decay to the latter, which 
will be reflected by the fact that the energy of dense matter on is higher than that on T^. 
As indicated in the introduction, the Skyrmion analysis shows that the energy on is only 
few percent that of T^. So we expect the current approximation to capture the essentials 
of dense matter in holographic QCD. Specifically, the nature and strength of the attraction 
and repulsion in dense matter. Indeed, this will be the case as we will detail below. 

Now, consider Nf probe D8-branes in the D4-branes background. With U{Nf) gauge 
field Am on the D8-branes, the effective action consists of the DBI action and the Chern- 
Simons action 



Sbs — Sbbi + Scs , 

Sbbi = -Ts J d^x e"^ tr V-det(5Miv + 27ra'FMiv) , (4) 

where Ts = l/((27r)^/3), the tension of the D8-brane, Fmn = Qm^n — On Am — i [Am, An] 
(M, = 0, 1, ■ • • ,8), and qmn is the induced metric on the D8-branes 

'^^^^ = [rJ ^'"^^^ + ^'"^^'^ + + Vu) ^''^^^ ' 

g^a = GrrdaTdaT + GuudaUdaU , (7) 



where Gmn refer to the background metric (1) and the profile of the D8 brane is parame- 
terized by U{a) and T{a). 

The gauge field Am has nine components, ^Iq, A — ^1,2,3, -^ai— ^4), ^-nd Aa{a — 
5, 6, 7, 8, the coordinates on the -S"^). We assume 

Ao = Ao{a) e U{1) , (8) 

(A = A = ^)) e su(Nf) , (9) 

A, = . (10) 
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Then the action becomes 5-dimensional: 



'DBi = °" tr / dte^da U 



s 



1/2 



'^cs = ^tr JaaFAF, (12) 

where €3 is the volume form of 5"^ space and the indices i, j, k{e {^p, 9, 0}) are raised by the 
metric g''^ defined by 



' J^^'D^^sin^V'' J^^sin^V'sin^ J ' ^^^^ 



3 Instanton in x 

Only Aq will be determined dynamically in the given instanton background Ai, A^. The exact 
background instanton solution is unknown. Thus we start with an approximate solution 
which is the SU(2) Yang-Mills instanton solution in the space with metric, 

ds^ = da^ + Oi^dQa . (14) 

This metric is different from our metric in (6) and (11), where there are warping factors. 
Furthermore our action is the nonlinear DBI action and not a Yang-Mills action. However 
it can be shown that the Yang-Mills instanton in the space (14) is the leading order solution 
of 1/ A expansion of the full metric and the DBI action as shown in [3] . So the solution can 
be used in the leading order calculation. 

We summarize here the (anti) self dual instanton solution obatained in [15]. Using the 
ansatz , 

A^f{a)U-^dU, C/ = cosV^ + ir„f"(^,0)sinV^ , (15) 
we get the field strength, in terms of vielbein whose relation to the coordinate ip, 9, 4>, is 
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specified in [15], 



e° A e" + - 



(16) 



where we used La — U ^daU — Ta/"^- If we require (anti) self-duality, 

2(/^ - /) 



3? 



(17) 



1 



(18) 



(19) 



then / is determined as 

1 _|_ g=F2(<7-o-o)/3l ' 

SO the field strength of one (anti) instanton solution is 

i^^ = (9./±)|(e°Ae«±^6V'Ae'^). 

4 D8 brane plus Instanton 

Now that we have the background instanton configurations, the remaining dynamical vari- 
ables are r and Aq. However it can be shown that (?o-r(a) = is always a solution of the 
Euler- Lagrange equation regardless of the gauge field. For simplicity we will work with this 
specific configuration so that the only dynamical variable is Aq. Let us parameterize Aq by 
Z defined as 



, K=l + Z^ . 



KK 



(20) 



Then the field strength is expressed in terms of Z and the dimensionless radius 3^ = ^R/C/kk, 



Fza = 77^/' 



ab 



Ukk :k 

2 J 



KK 



■31 



(21) 
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where /' = dzf- The instanton configuration is 



1 + g=F2(2-^o)/3^ 



(22) 



The DBI action reads 



S^^, = tr / dtesdZK'/' 

3g, J 



KK 



R 



+ ( ^] (2W)2F|„ + ^(27ra')^(e„,ei^„zF,e)' 



1/2 



(23) 



where /Iq = dzAo. We are using the same vielbein coordinates as (16). Since the instanton 
size (3?) is of order 0(A~^/^) we define a new dimensionless parameter IR, which is order of 
(AO), as 



31 



KK 



(24) 



Furthermore we rescale the coordinate and the instanton field strength for a systematic 1/ A 
expansion 



A-^/^x" , Z \-^I^Z , t^t , 

Fab AFab , FaZ ^ AF^^ , Aq ^ Aq , 

K^{l + Z')^fl + jZ''\^K,, 



so all coordinates and gauge fields become of order of 0(A°) and we 
By using the instanton solution (22) we get 



S- 



DBI 



1 + 



tr / dte-^dZK 
'KK J 
3^2 ,'2 



1/3 



-K 



1/3/ 1 3 TT ^1/3. 



4M|kC/|k ' A4Mkk 

2 



4/3 



1/2 



(25) 



(26) 
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If we let Ukk — -^kk ^"^^ simplicity, then the DBI action yields 



S- 



DBI 



-dNA I dte.dZ ^/ <; 1 + kI^'F^ - jKl^'iA^r ] {i^f + kI^'F^] 



—dNcX I dte^dZ 



1 



3A 2A 



{A,y + 0((1/A)^) 



, (27) 



where 



d = ^^^^ 



. _ , Jin = Tin , F^ = J 

3%5 ' " 2Mkk 4 



J = 



_ /'2 sech^(Z/Dl) 1 



5i 



dzX = dz- 



~4 ~3 

tanh(Z/^) ^1 + ^sech2(Z/^)^ 



60? 



.sgn(Z) 



(28) 



The Chern-Simons action does not change by the recaling (25) , and it is order of A*^. 
With the instanton solution (22) the Chern-Simons action reduces to 



N f N 
Scs = ;T7^tr AAF AF 



247r2 j 87r2 
= cN^ j dte^dZAoF'^ , 



tr / dte^dZAo^ (2AM^^^ 



where 



(29) 



4M^K 
3%^ 



(30) 



It also can be written as 



Scs = ^N^^ / dZAodzX ^ , for yio = 1 , 



(31) 



which confirms that the field configuration (21), and (22) describe the single (anti) self dual 
instanton since Scs corresponds to Nc x the Pontryagin index when yio = 1. 
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5 Equation of State in 1/A 

The equation of state of cold holographic matter is the energy following from the action 
functional. The total action up to order of A° is 



S = J dtesdZ {£jY)Bi + ^cs) 
= -dNc / dte^dZ 



+ cN^ / d^xdZAoF' . 



(32) 



where Aq is an auxiliary field with no time-dependence that can be eliminated by the equation 
of motion or Gauss law, 



(33) 



with 



n = ^ = dNcAo , 



OA, 







(34) 



The integral of the equation of motion with F in (28) is 



n(oo) 

where we have set 11(0) = 0. 
The energy of one cell is 



£cell — 



u{z) = n(oo) 



tanh(Z/:Ji) ^1 + ^sech2(Z/3?)^ 



^KK -'"c 

547r3 ^3 ' 



€3dZ{Lr)Bi + ^cs 
= dNr I €3 dZ 



_ 72 _ 
XF^ + —F^ 
3 



= dNc / eg dZ 



1 



AF^ + —F^ 



2 {dNc] 

1 



2 {dNc 



€3 dZAoW , 
£3 Ao{Z)U{Z) 



(35) 



(36) 
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Thus the energy density (e) of the crystaUine structure is 



iV£ 



cell ^cell 



dN^J dZ 



~2 ^'?2 1 n2 



-yio(^)n(z) 



(37) 



where A?" is the total number of baryons(cells) and V is the total volume which is approx- 
imated by N J €3. Interestingly the second term in (37) is equal to iib^b since the density 
and the baryon chemical potential is given by 



riB 



respectively and because 



liB = NcAo{oo) . 



(38) 



Ao{Z)U{Z) 



2yio(oo)n(oo) = A^c^o(oo; 



1 



= fJ'BriB 



(39) 



Notice that is of order of (A)° from (24) so the baryon density is of order (iVcA)°. Since 
the action is finite and concentrated in a finite size, we can restrict the integral to the region 
Z < Zc and expand the action in 1/A. 



e = dNr 



= Mo 



dZ 



~2 z^?=.2 1 n2 



(40) 
(41) 



where 



Mo = 87r^ft;MKK 
(tt^ - 6)M, 



a 



2 

KK 



36(27r2)2/3 



6 = 



2167r3 ' 
367r4 



2M3^ 



(42) 



A few remarks are in order. 



1. Zeis introduced as an arbitrary cut-off which is bigger than the instanton size. However 
in section 8, we will argue that Z^. should be identified as baryon size by explicitly 
contructing the Skyrmion out of the instanton Therefore it is not an arbitrary number. 



2. Even in the case the instanton size is small, the baryon size on the boundary is not. 
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It is of order {NcX)° and large in units of Mkk- This point is important. While the 
instanton size in bulk is of the order of the string length and thus small as 1/-\/A in 
units of MfCK, its image on the boundary is a large Skyrmion. 

3. The position of the instanton Zq in the conformal direction is set to zero by parity. 

The various density contributions in (41) can be understood from the zero density and 
finite instanton calculation discussed by Sakai and Sugimoyo to order NcX'^. For that, we 
recall that the energy balance for a holographic instanton with fiat directions reads 
schematically as [3] 

7Ve(^AAp2 + B^^ (43) 

leading to an instanton size in bulk of order p ~ (B/A)^/^/-\/A. The Coulomb repulsion B 

is 10^ times the gravitational attraction A resulting into a size that is of order p ~ 

This parametrically huge repulsion results in a stiffer equation of state in holographic QCD. 

The hnear term in ub in (41) is just the topological winding of the U{Nf) flavored 
instanton in D8 on 5"^ due to the self duality of the instanton configuration. It is leading and 
of order NcX. Geometry is unaffected by matter. A point-like instanton in bulk corresponds 
to a very large Skyrmion on the boundary. The term of order n^^^ is of order A^^^A^. It 
corresponds to the attraction due to gravity in bulk at finite size. Indeed, the energy of this 

~2 

term is of order Xp^ = 3? , as in (43) favoring smaller and smaller instanton. The energy 
per volume for this term is of order l/Ji. Since in matter the cell size is of the order of the 
interparticle distance the n^^ follows. The term of order is also of order NcX^. It 

~2 

stems from the Coulomb repulsion in bulk which is of order 1 / Xp^ — l/'X since the instanton 
is static in 4-space (space-plus-conformal) . This contribution is repulsive and favors larger 
size instanton. The corresponding energy per cell is of order [Zc/'X){l/'3l ), since the warping 
in the conformal direction is subleading in 1/A. The n\ contribution follows. 

For a Skyrmion with a size Zc <^ IR, (41) describes the low density regime. In this regime 
the use of the 5"^ x R instanton is hkely to give higher energy than a localized but flat 
instanton at the pole of say. Dilute holographic matter is made out of flat instantons 
with (41) providing an upper bound on the energy per unit volume. This phase breaks 
spontaneously chiral symmetry. In the point particle limit, the equation of state at low 
density was discussed in [8] 
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for low densities after re-scaling -s/Ane/A^/^ — > riB [8]. The point-like and flat space instanton 
contribution (44) at low density is lower in energy than (41) and therefore favored. This will 
be made more explicit below. The is absent in the point like limit (finite size effect). 

As the density is increased (or equivalently as ^ approaches down to Zc), there is a 
change in the equation of state (41). For Zc — Ji, 



with b changing to b' 



,« 1/3 , b' 5/3 



(45) 



3^(27r 



2\5/3 



The softening of the equation of state at higher density from to n^^^ follows from a 
transition from a dilute gas/liquid phase to a dense solid/crystal phase. This transition 
effectively restores chiral symmetry as we will show later. An estimate of the chiral transition 
density follows by comparing the n% term from (44) to the leading in (45) 

By setting €p — eg, the critical transition density follows 



4M 



n% = ^ ■ (48) 



6 Numbers 

To give some estimates of the numbers emerging from the current discussion, we first recall 
that in holographic QCD the mass of one baryon at next to leading order is not unique. We 
refer to [3] for a more thorough discussion. In particular, the baryon mass to order NcX'^ is 

Mb = Mo(i + ^) , (49) 

where c = 277r-\/2/15. Thus the interaction energy per unit volume for the dilute case is 

£;Diiute = ^ _ ^^^^ ^ ^ j^^^ V3 _ crzB + 6 n|) , (50) 
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Figure 1: The energy per unit volume: Zc = ^ (red) and = 3^ (blue). See text. 



while for the denser case it is , 

pDense _ _ „ ^^/r _ ^0 / 1/3 ^y!J>l'^\ 

Ant = £ — nBiViB — -T- yan^ — cub + o j . 



(51) 



For numerical estimates, we use Mkk = 500 MeV and Mq = 940 MeV for A''^ = 3 [3]. Our 
parameters are 

A ~ 53.2 , a ~ 0.095 fm"^ , b ~ 2172 fm^ , b' ~ 2039 W , c ~ 31 . 
The interaction energies are then 



^3 . ^Ujlute^QgYf^ 



^iDilute 
^int 

^3 . ^Dense 



lilt 

Dilute,? 
^int 



(GeVfm-=^) = 
(GeVfm-3) 



O.OOieSn}/^ - 0.548^5 + 192n| , 
0.00168nJ/^ - 0.548nB + 36.0 n^/^ 
-0.548ni3 + 60.3n| 



(52) 
(53) 

(54) 



where we used Zc = 5. Notice that due to the smallness of the coefficients of the first two 
terms, the last term is dominant even in the relatively small baryon density if it is not much 
smaller than 10~^fm~^. 



The results on S^: ^^t^"*" and E^^'"'"' are compared in Fig.(l) for = 5 ( = 5/M] 



KK 



Ifm with restored dimensions). The results on S"^ and R^ are compared in Fig. (2). In R^ 
the instantons are point-like or ~ oo [8]. The crossing from R^ to occurs at relatively 
small densities ~ 1.26 no with uq = 0.17 fm^^^'^ the nuclear matter density. 
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0.05 0.1 0.15 0.2 0.25 0.3 
Figure 2: The energy per unit volume: Zc = oo (red) and Zc = 'Ji (blue). See text. 



7 Equation of State in General 

The approximation of by suggested at the beginning of the paper was justified in 
way in the dilute limit or for small densities. Phenomenologically, we have found that the 
chiral phase transition from to x R occurs at few times nuclear matter density in 
holographic QCD, which is reasonable. The small size instantons dominate dense matter. 
This means that higher order corrections to both the DBI action and the starting D4 metric 
are important. While we do not know how to assess them, we now suggest that they may 
conspire to be small. Indeed, if we were not to expand the DBI plus CS actions, that is if 
we were to include only these class of higher order corrections our numerical results change 
only mildly. 

Consider, the total (DBI + CS) action is written as 



S = -dN^X J d^xdZ^A - jBiA'^y + cN^ J d^xdZJAo , (55) 



where 



A = K^/^ + /^-i/s J ^5/3 J T^V^ j2 (.^^ 

^ - + .^2 TTi ^\ ^ + tt2 ^\ ^ + i\/r2 ttG ' l^Oj 



3^TT^ ^ 3 . 2M|k 
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The equation of motion is 



n' = av, J , 



with 



n = 



The integral of of motion is 



u(z) = n(oo) 



tanh(Z/3^) ( 1 + hech^{Z/Ol) 



n(oo) = o - 

6^ 47r2([/KK3l)3 



VbU{oo) . 



The energy per cell is 

Sceii = —J e3dZ{Cjr)Bi + ^cs) 

I 



dNcX / 63 dZ 



B + 



- 63 dZAoW 



dN^X / ea dZ \ A 



AU^ 



XN^d^B 



- / es Ao{Z)U{Z) 



The energy density (e) of the crystalline structure is then 

NEcell £cell 



dNcX / dZ \ A + 



AU^ 

xmPB 



- Ao{Z)U{Z) 



where N is the total number of baryons (cells) and V is the total volume which 
mated hy N J 63. The second term in (37) is 



Ao{Z)U{Z) = 2yio(oo)n(oo) = A^c^o(oo) 



1 



= fJ'BriB 
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Figure 3: The energy per unit volume. 



where 



1 



/es 27r2(f/KK3i)3 27r2(v^D?)3 ' 
Since ub is fixed we may set \ib = 0. Then the energy density is 



(64) 



xmd^B 



K 



2/3 



(65) 



where we subtracted the vacuum value. Thus the interaction energy per unit volume is 

= e- ubMb . (66) 

In Fig. (3) we show the equation of state for the expanded and unexpanded actions. As ex- 
pected the corrections are of order 1/ A for a finite and small size instanton. The unexpanded 
energy is finite for any size of the instanton due to the gravitational warping factors which are 
subleading in 1/A after rescaling. The unexpanded results are similar to the expanded ones 
in the range of densities explored as it should. For extremely small hb-, -Emt ~ 0.00251n 



1/3 

B ' 



o 5/3 

however for reasonably large density ~ Ifm , E'int ~ 33.9 . This power is consistent 
and expected from the expansion in eq.(45). 

All general expressions in this section are consistent with the results quoted above to 
order 0(A°). Indeed, if for simplicity we set [/kk = ^kk with ^kk = 1, then A and B 
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reduce to 

12^2 

S = 6(1 + 36J) + 0(A-^) , (67) 



A^{1 + SbJf + Y^(l + 3&J)(2 + 36J) + 0(A-') , 
A o 



For example, by considering SbJ — and 11 = VbU, we can readily show that (65) reduces 
to (40) 



e ^ dN^ I dZ 



~2 1 n2 



3 2{dN,)\ 

8 Holographic Skyrmions from Instantons 

The S'^ X R instanton used in bulk has a very simple Skyrmion picture on the boundary. 
From (15) it follows that the gauge field at the boundary is A(oo, x) = U~^dU . Following [3] 
we note that U{x) is just the pion field at the boundary. When we have a cut-off in Z, we 
replace A[oo,x) by A{Zc,x). U is the boundary Skyrmion field originating from the bulk 
instanton. Thus U is just the holonomy of the hulk instanton along the conformal direction: 

U{x; Z,) = P exp[i ^z)] (68) 

Jo 

When the density is large and Zc ~ 3?, the instanton has a support covering the whole three 
sphere, therefore the resulting Skyrmion should be 

U{x) ~ (7(f) + iTaWix) = e^""''"(^''^)'^ , (69) 

which is the identity map as {ip, 6, (j)) are the canonical angles for the unit S'^ . The local 
Jacobian matrix for this map from to is J"* = dlV^/dxi = 1"*//?, proportional to the 
identity. The baryon density for this map is detJ/volS^ = 1/(27?^/?^) in agreement with bulk 
holography. The scalar field a{x) — cosip measures the chiral condensate and averages to 
zero on 



< qq >R3 



2 r 

< a{x) >53= - / dipsm'^ipcosip = . (70) 

Jo 



The xR instanton in (15) corresponds to a boundary Skyrmion on with restored chiral 
symmetry on the average. We should notice that the chiral condensation is p-wave over a 



16 



(a) 



3i 

Cb) 



Figure 4: Holographic Skyrmion on 5"^ on the boundary 

cell while the density in this case is approximately constant over a cell. But it is certainly 
not a constant. In fact this is a result consistent with ref. [7] where it was argued that there 
can not be an uniform distribution. In Fig. 4, we show schematically how a Skyrmion of size 
Zc looks on as a function of IR. (a) corresponds to the dilute phase with broken chiral 
symmetry, while (b) describes the dense phase with restored chiral symmetry. 

In previous section, Zc was introduced as a cut-off of the action bigger than the instanton 
size. Here we give interpretation of Zc as the size of the Skyrmion on the boundary. Note 
that the BPS instanton used in bulk in [3] for the description of a single baryon, yields 
a boundary Skyrmion as 

U{x)^Zclic-iT-xlic (71) 

with = Z^+x^+p^ and this is the analogue of the unit map (69) with tan'^ — x/^c- Notice 
that while the size of the instanton is p, the size of the Skyrmion is -^f Z"^ -\- p^. If p <^ Z^, 
Zc itself is the size of the Skyrmion, hence our interpretation above comes. Holography 
transmutes a small size instanton p in bulk to a large size Skyrmion on the boundary. 

At small densities with 0? ^ Zc, one can replace the spherical cell by a flat space and 
the map (71) is relevant, while at high density '5l< Zc the map (69) is relevant. On S"^ this 
is pictorially depicted in Fig. (4). Notice also that (a) has broken chiral symmetry while (b) 
has restored chiral symmetry effectively ( See eq. (70)) . Again, in this case, our x R 
instanton in bulk describes the high density phase in holographic QCD with restored chiral 
symmetry. At low densities the energy density is about n% as discussed by many [6] in 
qualitative agreement with our figure here. The term is sourced by Coulomb's repulsion 
in both cases. The description on carries larger energy density than on R'^ and is therefore 
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Figure 5: Holographic Skyrmion mass on S^: order A 



unfavorable energitically. It is favorable at higher densities. The transition occurs at about 
3? = Zc, or = l/(27r^ Z^), resulting into an energy density of n^^. The value of n% was 
estimated above. 

The determination of Z^ or equivalently the critical size of 51 depends on the energetics 
of the SS model. It is worth pointing that the single baryon mass analysis on S^' as discussed 
in [13] allows a considerable simplification of this issue when the Skyrmc model is used. Wc 
now note that this is justified in holographic QCD as small size instantons in bulk with 
p = "Ji/^/X map onto a large size Skyrmion on the boundary with Zc ^ p. So the small 
size instanton expansion in bulk maps onto the gradient expansion in 1/Zc on the boundary. 
Limiting the SS model on the boundary to the Skyrme model with fj^ and es fixed by 
holography yields the specifics of the Skyrmion on the boundary to order A. 

In Fig. (5) we show how the holgraphic Skyrmion mass on to order A changes with 31 
the radius of following [13]. The units of mass and length are respectively [3] 



= (AiV.)MKK^ (72) 



^ = (1/Mkk) v^A^ 

with b = 15.25 and L — 31. We note that the mass Mq — Stt^kMkk corresponds to the point 
0.95 at = oo which matches the unit map result as expected. In Fig. (6) we show the 
same curve to order 1/A. Here the energetics is determined in bulk as the chiral Lagrangian 
in the SS model is not known beyond the order A. Specifically, 
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Figure 6: Holographic Skyrmion mass on 5"^: order A° 



(73) 



with A = (tt^ - 6)/36 ~ 0.11 and B = {?>^ti'^)/^ ~ 1799. The units of mass and lengths are 



(B/A)V^ 
2\/AB 



~ 11.4 

- 27.8 



(74) 



The point 1.12 is the 1/A corrected mass (49) in these units. Finally, it is interesting to note 
that the holographic Skyrme model on yields naively the following equation of state 



e = Mo(nB + as Uji^ + hs n^^^) , 



(75) 



as first noted in the context of the canonical Skyrme model [16]. The rv^^ for the Skyrmion 
stems from the universal current algebra (VH)^ term which is attractive and scales as 1/p^ as 
opposed to 1/p from the finite size instanton in bulk. The for the Skyrmion stems from 
the repulsive Coulomb contribution per unit 3- volume (l/p)/p^ from the Skyrme term as 
opposed to the repulsive Coulomb contribution per unit 3- volume (1/p^)/ p'^ in the instanton 
in bulk. We recall that Coulomb's law in 1+D dimensions is 1/ p^~'^. 

At high density the naive scalings in (75) obtained at the boundary differs from (45) 
obtained in bulk in two essential ways: i) and hg are of order N^\^ on the boundary while 
their bulk contributions are of order N^/X; ii) the scaling with ub appears to differ by an 
extra (spatial) dimension, D = 3 on the boundary and D = 4 in bulk. These differences can 
be understood by noting that the size of the holographic Skyrmion is Zc- This means that 
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Figure 7: The energy per unit volume as a function of baryon density, for pure skyrmions, for 
the calculations of Bethe and Johnson [16, 17] and of Friedman and Pandharipande [16, 18], 
and for our instanton model based on Sakai-Sugimoto model for dilute and dense case. 



the chiral gradients Li = U~^diU are nearly zero on the boundary with f/ ~ 1, except on an 
the shell \x\ ^ Zc of thickness 1/VX to ensure that the topological baryon charge is finite ^. 
This renders as and bs in (75) effectively of order 1/A as noted in bulk. 

9 Comparison with Nuclear Models 

In Fig. (7) we compare the interaction energy (50) and (51) with other hadronic models 
including Skyrme's chiral model. Holographic matter is subtantially stiffer as explained 
through the energy budget in (43). The reason can be traced back to the fact that for a 
single baryon the repulsion already dwarfs the attraction in holographic QCD 
At high densities e in (51) is approximated as 

' - - 36nf (GeVfm-^) foriV. = 3, (76) 

and whatever Nf since the flavoured instanton in bulk is always 2x2. This behaviour is 
different from that of free massless quarks in D = 3 (£3) but similar to massive quarks in 
D = 3 (e'g). Specifically, 

£3 = ^^^!^^^^n^/^ ~ 5.52n^/^(GeVfm-3) for N, = 3, Nf = 2 , (77) 
^ K 33/571^/3 ^ ^ n'/'(GeVfm~=^) for A^, = 3, A^. = 2 . (78) 



^We note that for f7 ^ 1 the Skyrmion obeys the Faddeev-Bogomolnyi bound since the classical equations 
of motion are fulfilled. 
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So at strong coupling 




m 



KK 



Mkk 



88m 



iorNf = 2 , 



(79) 



independently of A and Nc- As chiral symmetry is restored in the high density phase, the 
comparison to the the free massive quark phase in D — 3 suggests that the mass m ~ Mkk/88 
is a chirally symmetric screening mass. While the chiral transition restores chiral symmetry 
it still confines baryons. 

10 Conclusions 

Wc have provided a holographic description of dense and cold hadronic matter using the 
brane model put forward by Sakai and Sugimoto [3]. At large Nc the matter crystallizes and 
can be treated in the Wigner-Seitz approximation on T^. For simplicity, the Wigner-Seitz 
cell was further approximated by in space leading to a simple instanton configuration on 
5*'^ X R with R the conformal space. The resulting equation of state at next to leading order 
in A shows a free quark behavior at high density, although the overall coefficient is cutoff 
sensitive and large resulting into a stiff euation of state. 

At high densities the gauge gradients are of order VX so the DBI action may not be 
enough to fix the brane dynamics at order NcX^ [3]. Also our simphfication of by 
while justified at low density, involves curvature corrections at high densities. However, we 
believe that the essentials of dense matter in holographic QCD are already exposed on 
with a small attraction leading and a large Coulomb repulsion leading n^^, where 5/3 is 
the power of non-relativistic fermion. It is interesting to notice that the coulomb interaction 
in the bulk counts the fermi statistics in the boundary. The repulsion is 10^ times the 
attraction resulting into a very stiff equation of state. Changing to will not affect the 
outcome quantitativaly we believe. Indeed, this is the case for dense Skyrmions [13]. 

The present work expands on the original ideas developed in [6]. Our calculations with 
finite size instantons are closer to those presented in reference [7] where finite size and 
homogeneous instantons were used through a variational estimate in R^xR. Their arguments 

1/3 

yield hb instead of the we have reported in the equation of state at next to leading order 
with our X R instanton. 

The inhomegeneous S^' x R description of the crystal suggests that at high density, chiral 
symmetry is restored on the average. Indeed, since the dual of the instanton cell is the 
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Skyrmion cell with a pion field restricted to in space. High density matter corresponds to 
small size where the pion field becomes just the unit map [13]. The corresponding chiral 
condensate on is seen to vanish as half of carries positive chiral condensate, while the 
other half carries negative chiral condensate so that on the average the chiral condensate 
is zero. This restoration of chiral symmetry is due to the formation of the crystal in the 
spatial direction in holographic QCD even though the D8-D8 configuration is still attached. 
In other words, the left and right D8 branes cease to talk to each other through the spatial 
directions not the conformal direction when they crystallize at large Nc- 

The present crystal analysis is classical in bulk. A quantum analysis including vibrational 
and rotational motion is needed. These corrections are subleading in 1/Nc and should be 
estimated for a more thorough phenomenological discussion. Also, the inhomogeneous phase 
can be probed approximatly by a dilute gas of instantons on allowing for a lower energetics 
than on S^. These issues and others will be discussed elsewhere. 
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